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Abstract 

The thermodynamic properties of hot nuclei are described within the canonical and microcanon- 
ical ensemble approaches. These approaches are derived based on the solutions of the BCS and 
self-consistent quasiparticle random-phase approximation at zero temperature embedded into the 
canonical and microcanonical ensembles. The obtained results agree well with the recent data 
extracted from experimental level densities by Oslo group for 9 Mo, 98 Mo, 162 Dy and 172 Yb nuclei. 
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I. INTRODUCTION 



Thermodynamic properties of highly excited (hot) nuclei have been a topic of much in- 
terest in nuclear physics. From the theoretical point of view, thermodynamic properties of 
any systems can be studied by using three principal statistical ensembles, namely the grand 
canonical ensemble (GCE), canonical ensemble (CE) and microcanonical ensemble (MCE). 
The GCE is an ensemble of identical systems in thermal equilibrium, which exchange their 
energies and particle numbers with the external heat bath. In the CE, the systems ex- 
change only their energies, whereas their particle numbers are kept to be the same for all 
systems. The MCE describes thermally isolated systems with fixed energies and particle 
numbers. For convenience, the GCE is often used in most of theoretical approaches, e.g. 
the conventional finite-temperature BCS (FTBCS) theory [1], and/or finite-temperature 
Hartree-Fock-Bogoliubov theory j^j]. These theories, however, fail to describe thermody- 
namic properties of finite small systems such as atomic nuclei or ultra-small metallic grains. 
The reason is that the FTBCS neglects the quantal and thermal fluctuations, which have 
been shown to be very important in finite systems 3-8|. These fluctuations smooth out 
the superfluid-normal (SN) phase transition, which is a typical feature of infinite systems as 
predicted by the FTBCS theory. 

Because an atomic nucleus is a system with the fixed particle number, the particle- 
number fluctuations are obviously not allowed. The use of the GCE in nuclear systems is 
therefore an approximation, which is good so long as the effect caused by particle-number 

'he CE and MCE are often used in extending the exact solutions 



fluctuations are negligible 
of the pairing Hamiltonian 8MlO| to finite temperature, whereas the CE is preferred in the 
quantum Monte-Carlo calculations at finite temperature (FTQMC) [l]} Q]- However, it 
is impracticable to find all the exact eigenvalues of the pairing Hamiltonian to construct 
the exact partition functions for large systems. For instance, in the half-filled doubly-folded 
multilevel model (also called the Richardson model) with N = Q with Q being the number 
of single-particle levels and N - the number of particles, this cannot be done already for 
N > 14 [8|, [9]. Meanwhile, the FTQMC is quite time consuming and cannot be applied to 
heavy nuclei unless a limited configuration space is picked up. It is worth mentioning that 
the pairing Hamiltonian can also be solved exactly by using the Richardson's method, i.e. 
solving the Richardson equations. Using this method, the lowest eigenvalues of the pairing 



Hamiltonian can be obtained even for very large systems, e.g. with N = Q = 1000 (See, 
e.g., Ref. Q). Nonetheless, these lowest eigenstates (obtained after solving the Richardson 
equations) are not sufficient for the construction of the exact partition function at finite 
temperature since the latter should contain all the excited states, not only the lowest ones. 
In principle, the CE-based approaches can also be derived from the exact particle-number 
projection (PNP) at finite temperature on top of the GCE ones However, this method 
is rather complicated for application to realistic nuclei. 



The static path plus random 



P 



rase approximation (SPA + RPA) with the exact num- 



ber parity projection CSPA(p) 15j and the latter extension of the number projected SPA 



(NPSPA) 16j offer quite good agreement with the exact CE of the Richardson model as well 



as the empirical heat capacities of heavy nuclei. However, Ref. |15[ makes no comparison 



with experimental data, whereas Ref. 16( uses a thermal pairing gap, which is obtained from 
a direct extension of the odd-even mass difference to finite temperature. As has been pointed 
in Ref. js| such simple extension fails in the region of intermediate and high temperatures. 
In principle, the SPA can also be used to evaluate the MCE quantities based on the GCE 



ones by fixing the energy and particle number of the system [17J. However this method is 
still quite complicated for practical applications to realistic nuclei, especially the heavy ones. 
From the experimental point of view, the CE and MCE are usually used to extract vari- 
ous thermodynamic quantities of nuclear systems. This is carried out by using the nuclear 
level density, which can be experimentally measured at low excitation energy E* < 10 MeV. 
Within the CE, the measured level densities are first extrapolated to high E* using the back- 
shifted Fermi-gas model (BSFG). The CE partition function is then constructed making use 
of the Laplace transformation of the level density. Knowing the partition function, one can 
calculate all the thermodynamic quantities within the CE such as the free energy, total 
energy, heat capacity and entropy. The thermodynamic quantities of the systems obtained 
within the MCE are calculated via the Boltzmann's definition of entropy. Although several 
experimental data for nuclear thermodynamic quantities extracted in this way by the Oslo 



group have recently been reported [18l-l21|. most of present theoretical approaches, derived 
within the GCE, cannot describe well these data, which are extracted within the CE and 
MCE. Recently we have proposed a method, which has allowed us to construct theoreti- 
cal approaches within the CE and MCE to describe rather well thermodynamic properties 



of atomic nuclei 



221 ] . The proposed approaches are derived by solving the BCS and self- 



consistent quasiparticle RPA (SCQRPA) equations with the Lipkin-Nogami (LN) PNP for 
each total seniority S (number of unpaired particles at zero temperature) 23| . The obtained 
results are then embedded into the CE and MCE. Within the CE, the resulting approaches 
are called the CE-LNBCS and CE-LNSCQRPA, whereas they are called the MCE-LNBCS 
and MCE-LNSCQRPA within the MCE. The results obtained within these approaches are 
found in quite good agreement with not only the exact solutions of the Richardson model 
but also the experimentally extracted data for 56 Fe isotope. The merit of these approaches 
reside in their simplicity and feasibility in the application even to heavy nuclei, where the 
exact solution is impracticable and the FTQMC is time consuming. The goal of present 
article is to apply the above-mentioned approaches to microscopically describe the recently 
extracted thermodynamic quantities of 94,96 Mo, 162 Dy and 172 Yb nuclei. 

The article is organized as follows. The pairing Hamiltonian and the derivations of the 
GCE-BCS, CE(MCE)-LNBCS and CE(MCE)-LNSCQRPA are presented in Sec. HI The 
numerical results are analyzed and discussed in Sec. 1111} whereas the conclusions are drawn 
in the last section. 



II. FORMALISM 

A. Pairing Hamiltonian 

The present article considers the pairing Hamiltonian 

H= y j e k al a a k(T - G) j al + al_a k ^a k/+ , (1) 

fctr=± kk' 

where a\ a and a ka are particle creation and destruction operators on the kth. orbitals, respec- 
tively. The subscripts k here imply the single-particle states in deformed basis. This Hamil- 
tonian describes a system of N particles (protons or neutrons) interacting via a monopole 
pairing force with constant parameter G. The pairing Hamiltonian (JT|) can be diagonalized 
exactly by using the SU(2) algebra of angular momentum [10] . At finite temperature T ^ 0, 
the exact diagonalization is done for all total seniority or number of unpaired particles S 
because all excited states should be included in the exact partition function. Here S = 
0, 2, ... N for even-A systems, and 5 = 1,3,... N-l for odd-A systems. For a system 
of A particles moving in fl degenerated single-particle levels, the number riExact of exact 
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eigenstates £^ xact (is =1, . . . , ^Exact) obtained within exact diagonalization is given as 

n Exa ct = C s x C n~ S a . (2) 

' * iV pair o 

5 

which combinatorially increases with AT, where = m!/[n!(m — n)!] and iV pa i r = iV/2 jg]. 
Therefore, the exact solution at T 7^ is impossible for large N systems, e.g. N > 14 for 
the half-filled case (N = Q), because of the huge size of the matrix to be diagonalized. 



B. GCE-BCS 



The well-known finite-temperature BCS (FTBCS) approach to the pairing Hamiltonian 
(CQ) is derived based on the variational procedure, which minimizes the grand potential 

Q = (H) — TS — XN so that 5Vl = , (3) 

where S is the entropy of the system at temperature T. The chemical potential A is a La- 
grangian multiplier, which can be obtained from the equation that maintain the expectation 
value of the particle-number operator to be equal to the particle number N. The expectation 
value (O) denotes the GCE average of the operator O [6j (the Boltzmann's constant k B is 
set to 1), 

Tl \ 0e -(3(H-\N)1 x 

<0> 5 p = f- (4) 

The conventional FTBCS equations for the pairing gap A and particle number N are then 
given as 

A = G^2u k v k (l-2n k ) , N = 2^2[(l-2n k )vl + n k ] , (5) 

k k 

where the Bogoliubov's coefficients u k , v k , the quasiparticle energy E k and the quasiparticle 
occupation number n k have the usual form as 



2 V E k J 2 \ E, 

1 



E k = ^(e k -Gvl-\y + A\ n k = YV -^. (6) 

The systems of Eqs. (jSJ) and (jSJ) are called the GCE-BCS equations. The total energy, heat 
capacity and entropy obtained within the GCE-BCS are given as 

£ = 2 [(! - 2n k )v 2 k + n k ] - ^ - G J> - 2n*)^ , 



C = — , S = -2 ^2 ["felnrife + (1 - ra fc )ln(l - n fc )] 

k 



(7) 



C. CE-LNBCS 



Different from the GCE-BCS, the CE-LNBCS is derived based on the solutions of the BCS 



equations combined with the Lipkin-Nogami particle-number projection (PNP) 24| at T = 
for each total seniority S of the system. When the pairs are broken, the unpaired particles 
denoted with the quantum numbers ks block the single-particle levels k. As the result, 
these blocked single-particle levels do not contribute to the pairing correlation. Therefore, 
the Lipkin-Nogami BCS (LNBCS) equations at T = can be derived by excluding these ks 
blocked levels. These equations are given as 

A LNBCS (M = G £ u k v k , N = 2j2v 2 k + S, (8) 

kj^k s k^k a 



where 



1 A , e k - Gvl - X(k s ) \ , If e k - Gvj - \{k s ) , 

Vk+k s - 2 \1 + ^ ) » %* fl - 2 ^ ^ 1 ' (J) 



S Mfes = ^/[c* - Gt£ - A(Afr)]' + [Alnbcs (a , 5 ) ]2 ) (10) 

w, s x /, s ^ /, \/»r ,\ * /. \ G^k^k s U k V k^ k '^ k ' U k'vl> — J2k^k s U k V k 

X(k s ) = + 2\ 2 (k s )(N + 1) , A 2 (A; 5 ) = j^j^ o% 5 2 v -44 • 

* {l^kj^ks U k u k) 2^ik+k s U k u k 

(11) 

As for the blocked single-particle levels, k = ks, their occupation numbers are always equal 
to 1/2. Solving the systems of Eqs. ([8]) - (1111) . one obtains the pairing gap 
quasiparticle energies E k and Bogoliubov's coefficients u k and v k , which correspond to each 
position of unpaired particles on blocked levels ks at each value of the total seniority S. 
There are tt-lnbcs = ^2s configurations of ks levels distributed amongst fl single-particle 
levels at each value of S, which is also the number of eigenstates obtained within the LNBCS. 
The LNBCS energy (eigenvalue) for each configuration is then given as 

[ALNBCSfi, \12 _ 

£Lnbcs = 2 J2 e k vl + J> s - 1 G { S)1 -Gj^vt- AX 2 (k s ) £ u\v\ . (12) 

kjtk s k s k^k s k^k s 

The partition function of the so-called CE-LNBCS approach is constructed by using the 



LNBCS eigenvalues £™ BCS as [22 1 



"LNBCS 



ws(/o = I> £ ^ > ( 13 ) 

S is=l 



where ds = 2 s is the degeneracy. Knowing the partition function (!T3|) . we can calculate all 
thermodynamic quantities of the system such as the free energy F, entropy S, total energy 
£, and heat capacity C as follows 



BF BP 
F = -T\nZ(T), S = -—, £ = F + TS, C = —. (14) 

The pairing gap is obtained by averaging the seniority-dependent gaps A^ s NBCS = A LNB (ks) 
at T = in the CE by means of the CE-LNBCS partition function (1131) . namely 

"LNBCS 

A - 1 A LNBCS — /3fi^ BCS f.K\ 

^CE-LNBCS — ~ Z, a S 7, A i s e S • \ lti ) 

^LNBCS a 



D. CE-LNSCQRPA 

As mentioned previously in sec. Ill a complete CE partition function should include 
all eigenstates. The LNBCS theory (at T = 0) produces only the lowest eigenstates. For 
instance, for even (odd) N there is only one state at S = 0, which is the ground state. 
For S > there are also excited states in even (odd) systems, whose total number nLNBCS 
is much smaller than nExact- Consequently, the results obtained within the CE-LNBCS 
can be compared with the exact ones only at low T because at high T, higher eigenstates 
(excited states), which the LNBCS theory cannot reproduce, should be included in the 
CE partition function. This can be done by going beyond the quasiparticle mean field 
by introducing the SCQRPA with Lipkin-Nogami PNP (LNSCQRPA), which incorporates 



not only t 



QRPAJ221 
Refs. 
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re ground states but also the pairing vibrational excited states predicted by the 
he derivation of the LNSCQRPA equations has been presented in details in 



23 



251 ]. so we do not repeat it here. The LNSCQRPA formalism at T = for 
each total seniority S is proceeded in the same way as that of the LNBCS described in the 
previous section, namely the LNSCQRPA equations are derived only for the unblocked levels 
k ks, whereas the levels, blocked by the unpaired particles k = ks, do not contribute to 



the pairing Hamiltonian. The SCQRPA equations at T = has been derived in Ref. 23], 
whose final form reads 



A B \ I XI \ [XI 




B A I \ l \ -Y£ 



uul , (16) 



The SCQRPA submatrices are given as 
A kk > = 2 



b k + 2q k k> + 2 toK 1 - v k") - p- ( S 4fc»(0|^I»A|0) 
fc" fc ^ fc" 

- 2 V ^ fc »(0|A»A10))1 <W + d kk ,y/V k V k , + 8g fcfc , (5| *^ft , (17) 
77 /J V^k^k' 

h kk > + dkk " (01 A» A |0) + 2 h kk" (0\Al„ A k |0) 



fc 

-Bfcfc' = — 2 



+ 2^v / ^ + 8gfcfc^ (0| ^^° ) , (18) 
where d kk », h kk », and g^fc' (ah k ^ &£•) are functions of u k , v k , e k , A and G as given in Eqs. 



(13), (15), (17) and (18) of Ref. |23|. The screening factors (0|^A'|0) and (0\A k A k >\0) 
with A^ = or k a_ k being the creation operator of two-quasiparticle pair are given in terms of 
the SCQRPA amplitudes A£ and y v k as 

(oi4A'io) = vWEw , (o\A k A k ,\o) = VtPkfiPdY,**?* ■ ( 19 ) 

V V 

where (0| . . . |0) denotes the expectation value in the SCQRPA ground state. The ground- 
state correlation factor D k is expressed in term of the backward-going amplitudes y k as 
V k = [1 + 2 'YlvO^kYY 1 "with the sum running over all the SCQRPA solutions v. 

After solving the LNSCQRPA equations (jHJ), ( |T6i) - ( fl8l) for each total seniority S, we 
obtain a set of eigenstates, which consists of C§ lowest eigenstates (the ground state at 
S =0 or 1), as well as higher eigenstates (excited states) on top of these lowest ones, which 
come from the solutions of the LNSCQRPA equations, whose eigenvalues are ool S ^ [y = 
1, . . . Q— S^. As the result, the total number of eigenstates obtained within the LNSCQRPA 
is given as 

^LNSCQRPA = £)Cgx(fi-S). (20) 

s 

Consequently, the so-called CE-LNSCQRPA partition function is calculated as 

"LNSCQRPA 

Zlnscqrpa(/3) = e 13 > ( 21 ) 

S i 5 =l 



1 The SCQRPA has altogether O— S+l solutions with positive energies. However the lowest one corresponds 
to the spurious mode, whose energy is zero within the QRPA. Therefore it is excluded in the numerical 
calculations. 
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which is formally identical to the CE-LNBCS partition function (fT3|) . but the LNBCS eigen- 
values £j^ BCS are now replaced with cqrp ^ p rom partition function, the thermo- 
dynamic quantities obtained within the CE-LNSCQRPA are calculated in the same way as 
those in Eq. ( 1T4~|) . Although the number ^lnscqrpa of the LNSCQRPA eigenstates is larger 
than nLNBCS; it is still much smaller than nE xact . This most important feature of the present 
method tremendously reduces the computing time in numerical calculations for heavy nuclei. 
As an example, we show in Tabled the number of eigenstates and the total executing time 
(the elapsed real time) for the exact diagonalization of the pairing Hamiltonian, CE-LNBCS 
and CE-LNSCQRPA calculations within the Richardson model at several values N of par- 
ticle number, which is taken to be equal to number Q of single-particle levels (the half-filled 
case). This table shows that the execution time within the LNSCQRPA (LNBCS) is shorter 
than that consumed by the exact diagonalization by about two (four) orders. 



TABLE I. Number of eigenstates and computation time for the exact diagonalization of the pairing 
Hamiltonian as well as the numerical calculations within the CE-LNBCS and CE-LNSCQRPA for 
the doubly-folded equidistant multilevel pairing model at several values of N = f2. The computation 
time is estimated based on a shared large memory computer Altix 450 with 512GB memory of 
RIKEN Integrated Cluster of Clusters (RICC) system. 





Number of eigenstates 


Computation time 


N 


Exact 


LNBCS 


LNSCQRPA 


Exact 


LNBCS 


LNSCQRPA 


10 


8953 


512 


2560 


1 hr 


1 sec. 


10 sec. 


12 


73789 


2048 


12288 


10 hrs 


10 sec. 


1 min. 


14 


616227 


8192 


57344 


24 hrs 


1 min. 


10 min. 


16 


5196627 


32768 


262144 




10 min. 


1 hr 


18 


44152809 


131072 


1179648 




1 hr 


3 hrs 


20 


377379369 


524288 


5242880 




3 hrs 


10 hrs 
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E. MCE-LNBCS, MCE-LNSCQRPA 



The MCE entropy is calculated by using the Boltzmann's definition 



S(£) = \nW(£) , W(£) = p{£)5£ 



(22) 



where p(£) is the density of states. In the LNBCS (LNSCQRPA) approaches, W(£) is 
the number of LNBCS (LNSCQRPA) eigenstates within the energy interval (£,£ + 5£) js]. 
Knowing the MCE entropy, one can calculate the MCE temperature as the first derivative 
of the MCE entropy with respect to the excitation energy £, namely 

(23) 



T 



dS{£) 



d£ 



The corresponding approaches, which embed the LNBCS and LNSCQRPA eigenvalues into 
the MCE, are called the MCE-LNBCS and MCE-LNSCQRPA, respectively 



F. Level density 



The inverse relation of Eq. (122"]) reads 



p(£) = e s{£) /5£ 



(24) 



which can be used to calculate the density of states p{£ ) from the fitted MCE entropy. 
Within the CE, the density of states p{£ ) is calculated by using the method of steepest 



descent to find the minimum of the Laplace transform of the partition function 26]. As 
a result the density of states p(£) at temperature T = /3q , which corresponds to this 
minimum, is approximated as 



p{£) w Z%)e^ £ 



2tt 



d 2 lnZ((3 



-1/2 



2vr 



d£ 
Wo 



-1/2 



(25) 



where Z({3o), S(£) and £ are the CE partition function, entropy and total excitation energy 
of the systems, respectively. The density of states p(£) is obtained within the CE-LNBCS 
and CE-LNSCQRPA by replacing the partition function Z in Eq. (|25p with that obtained 
within the CE-LNBCS in Eq. flTJ and CE-LNSCQRPA in Eq. (J2U). 

At finite angular momentum J, in principle, the approach of LNSCQRPA plus angular 



momentum, which has been proposed by us in Ref. 



27| . should be used to calculate the 



10 



angular-momentum dependent level density p(£, M) with M being the ^-projection of the 
total angular momentum. In this case the former doubly-degenerated quasiparticle levels 
are resolved under the constraint M = Yl,k m k( n k ~ n k) with the quasiparticle occupation 
numbers rit, which are described by the Fermi-Dirac distribution np FD = {exp[(3(Ek =p 
7?7ifc)] + within the non-interacting quasiparticle approximation, where is the spin- 
projections of the kth single-particle state \k, im^), Ek is the quasiparticle energy, and 7 is 
the rotation frequency. Knowing p{£ , M), one can find p(S, J) = p(8, M = J) — p(£, M = 
J + 1) in the general case, where the total angular momentum J is not aligned with the 



z- axis 



281 ] . The total level density ptot(£) and experimentally observed level density p obs (£), 



are then defined as 



29] 



p tot (£) = ^(2J + l)p(£, J) , Pobs (£) = ^2p(£,J) ■ (26) 
j j 

The empirical entropy S b s (£) is extracted from the observed level density p bs(£) in the 

same way as in Eq. ( |22i) . replacing p(£) with p b s (£), namely 

S obs (£) = HPobs(S)5£] , (27) 

Because the present article considers non-rotating nuclei at low angular momentum, we 
assume that p(£, J) ~ p(£,0) = p(£). Therefore, by fitting the MCE entropy S(£) in Eq. 
(|22|) to the experimentally observed entropy S b s (£) in Eq. (1271) . i.e. S(£) ~ S b s (£), and 
inverting the obtained result by using Eq. (12^1) . what we get is actually the level density 
comparable to the experimentally observed one, p bs{£) = exp[S(£)]/8£. This means that 
the density of states p(£) calculated by using Eq. ( 1241) or Eq. ( |25l) without taking into 
account the effect of finite angular momentum is identical to the level density like p b s (£), 
not the total level density Ptot(£), because of the absence of the factor (2J + 1). 

III. ANALYSIS OF NUMERICAL RESULTS 

The proposed approaches are used to calculate the pairing gap, total energy, entropy and 
heat capacity within the CE and MCE for a number of heavy isotopes, namely 94,98 Mo, 162 Dy 
and 172 Yb @. The single-particle energies are taken from the axially deformed Woods-Saxon 



See, e.g. Fig. 1 of Ref. 



22| and Appendix A of the present article for the accuracy of the present 



approaches in comparison with the exact solutions of the Richardson model. 
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potential with the depth of the central potential 



3Q| 



V = V n 



1 ± k- 



N + Z 



(28) 



where Vq = 51.0 MeV, k = 0.86, whereas the plus and minus signs stand for proton (Z) and 
neutron (N), respectively. The radius tq, diffuseness a, and spin-orbit strength A are chosen 
to be ro = 1.27 fm, a = 0.67 fm and A = 35.0. The quadrupole deformation parameters 02 
are estimated from the experimental B(E2; 2^ — > 0i), which are 0.15, 0.17, 0.281 and 0.296 



for 94 Mo, 98 Mo, 162 Dy and 172 Yb, respectively 21|. The pairing interaction parameters G are 



adjusted so that the pairing gaps for neutrons and protons obtained within the LNSCQRPA 
at T = and S = reproduce the values extracted from the experimental odd-even mass 
differences, namely ~ 1.2, 1.0, 0.8 and 0.8 MeV for neutrons, and Az — 1.4, 1.3, 0.9 
and 0.9 MeV for protons in 94 Mo, 98 Mo, 162 Dy and 172 Yb, respectively. 

It is well-known that pairing is significant only for the levels around the Fermi energy. 



Therefore, within the CE, we apply the same prescription proposed in Ref. [12| to calculate 
the CE partition function for medium and heavy isotopes. According to this prescription, we 
calculate the LNBCS and LNSCQRPA pairing gaps in the space spanned by 22 degenerated 
(proton or neutron) single-particle levels above the doubly- magic 48 Ca core for Mo isotopes, 
whereas the same is done on top of the doubly- magic 132 Sn core for Dy and Yb nuclei. The 
obtained partition function is then combined with those obtained within the independent- 
particle model (IPM) by using Eq. (15) of Ref. 12J, namely 



\nZl = \nZ' u tr + \nZ' sp - \nZ' sp tr , (29) 

where Z' utr = Z Uitr e l3£o is the excitation partition function with respect to the ground state 
energy So with Z Vitr being the CE partition function obtained within the LNBCS [Eq. ffT3l) ] 
or LNSCQRPA [Eq. (1211) ] for 22 degenerated single-particle levels around the Fermi energy. 
Z' is the CE partition function obtained within the IPM [See e.g. Eq. (8) of Ref. 12]] 
for the space spanned by the levels from the bottom to iV = 126 closed shell, whereas Z' sp tr 
is the same partition function but for the truncated space spanned by 22 levels around the 
Fermi energy. 
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FIG. 1. (Color online) Pairing gaps A and heat capacities C obtained within the CE as functions 
of T and entropies S obtained within the MCE as functions of E* for 94 Mo (left panels) and 98 Mo 
(right panels). In (a) and (d), the solid and dash-dotted lines denote the pairing gaps for protons 
and neutrons, respectively, whereas the thin and thick lines correspond to the CE-LNBCS and CE- 
LNSCQRPA results, respectively. In (b) and (e), the thin and thick lines stand for the CE-LNBCS 
and CE-LNSCQRPA results, whereas the thin and thick dash-dotted lines depict the experimental 

rl 



results taken from Refs. [20] and 



21], respectively. Shown in (c) and (f) are the MCE entropies 



obtained within the MCE-LNBCS (rectangles), MCE-LNSCQRPA (triangles), and extracted from 



experimental data (circles with error bars) of Ref. 



A. Results for molybdenum 



Shown in Fig. [T] are the pairing gaps, heat capacities and entropies for 94 Mo [Figs. Q] 
(a)-HI(c)] and 98 Mo [Figs. Q] (d)-Q] (f)] obtained within the CE(MCE)-LNBCS and CE(MCE)- 
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LNSCQRPA versus the experimental data from Refs. [20JJ and [211 ] . There is a clear dis- 
crepancy in the heat capacities extracted from the same measured level density in these two 
papers [Figs. d](b) and[T](e)]. The heat capacity, extracted in Ref. [21], clearly shows a pro- 
nounced peak at T ~ 0.7 MeV for both 94 Mo and 98 Mo, whereas the corresponding quantity, 
extracted in Ref. 20|, shows no trace of any peak. The source of the discrepancy comes 



from the difference in the scale of excitation energy E*, which was used for extrapolating 
the measured level density before evaluating; the CE partition function using the Laplace 
transformation of the level density. In Ref. 20J , the level density is extrapolated up to E* ~ 



40 - 50 MeV, whereas in Ref. [2l| this is done up to E* = 180 MeV. Given that all the 
excited states should be included in the partition function, the energy E* ~ 40 - 50 MeV 



used in Ref. 



20[ seems to be too low, which might affect the resulting heat capacity. As 



Figs. CD (b) and [T] (e) show, the heat capacities predicted by the CE-LNSCQRPA are much 



closed to those obtained in Ref. 
tions for other nuclei 
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21] . They are also consistent with the FTQMC calcula- 



12| . It is important to emphasize here that quantal and thermal 



fluctuations within the CE-LNBCS (LNSCQRPA) indeed smooth out the SN phase transi- 
tion. As the result, the pairing gaps [Figs. [T] (a) and (d)] obtained for protons (solid lines) 
and neutrons (dash-dotted lines) within both CE-LNBCS (thin lines) and CE-LNSCQRPA 
(thick lines) do not collapse at the critical temperature T = T c of the SN phase transition, 
as predicted by the GCE-BCS, but monotonously decrease with increasing T. The neutron 
gap in Fig. □ (a) obtained within the CE-LNSCQRPA for 94 Mo (thick dash-dotted lines) 
is close to the three-point gap (dashed lines) obtained in Ref. [2lJ by simply extrapolat- 
ing the odd-even mass formula to finite temperature. As has been pointed out in Ref. js] 
such naive extrapolation contains the admixture with the contribution from uncorrelated 
single-particle configurations, which do not contribute to the pairing correlation. Therefore, 
to avoid obviously wrong results at high T, such contribution should be removed from the 
total energy of the system. Nonetheless, in the low temperature region (T < 1.3 MeV) as 
that considered here, where the contribution of uncorrelated single-particle configurations is 
expected to be small, the simple extension of the three-point odd-even mass formula to T ^ 
can still serve as a useful indicator. 



As has been discussed in Ref. [22|, at low E* the genuine thermodynamic observable 
is the MCE entropy because it is calculated directly from the observable level density by 
using the Boltzmann's definition (122]) . The experimental MCE entropies for 94,98 Mo are 
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FIG. 2. (Color online) Microcanonical entropy as function of E* obtained within the MCE- 
LNSCQRPA for 94 Mo using various values of energy interval 56. 

plotted in Figs. Q] (c) and ffl (f ) along with the predictions by the MCE-LNBCS and MCE- 
LNSCQRPA. These figures show that the MCE-LNSCQRPA results fit the available exper- 
imental data remarkably well. It is worth mentioning that the results obtained within the 
MCE-LNBCS (LNSCQRPA) are sensitive to the choice of energy interval 88 , which is used 
to calculate the number of accessible states W(£) in Eq. (J22l) . Figure [2] shows the entropies 
obtained within the CE-LNSCQRPA for 94 Mo using several values of 88 ranging from 0.2 
MeV to 1.0 MeV. It is clear to see from this Fig. [2] that the MCE entropies increase with 
increasing 68. In this respect, we found that the values of 88 = 1 MeV for 94 Mo and 0.7 
MeV for 98 Mo are reasonable to fit the experimental data. The reason for choosing large 
values of 88 for these two nuclei comes from the deficiency of the CE-LNSCQRPA(LNBCS), 
which includes only low-lying excited states. 

B. Results for dysprosium and ytterbium 

The results obtained for 162 Dy and 172 Yb are shown in Fig. |3j Similar to the results 
for 94 ' 98 Mo, the CE heat capacities and MCE entropies obtained within the CE(MCE)- 
LNSCQRPA for both 162 Dy and 172 Yb are in good agreement with the experimental data. 
The neutron and proton gaps obtained within the CE-LNBCS (LNSCQRPA) do not collapse 
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FIG. 3. (Color online) (a), (b), (e) and (f): Pairing gaps A, heat capacities C as functions of T 
obtained within the CE; (c), (d), (g) and (h): Entropies S and temperatures T as functions of E* 
obtained within the MCE for 162 Dy (left panels) and 172 Yb (right panels). Notations are the same 



as those in Fig. [TJ Experimental data are taken from Ref. 



at T = T c but decrease with increasing T and keep finite at high T even for the two heavy 
nuclei considered here. The peak in the experimental heat capacity near T = 0.4 MeV is 
seen in 172 Yb, whereas it disappears in 162 Dy. This is again due to the fact that the measured 
level densities for these two nuclei are extrapolated only up to E* =40 MeV instead of 180 
MeV as was done in Ref. [21 ] for other nuclei. This is confirmed by the heat capacities 
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obtained within the CE-LNSCQRPA (thick solid lines), which clearly show a peak around 
T = 0.4 MeV. 

In Figs. |3] (d) and |3] (h), one can see that the MCE temperatures, extracted from the 
experimental data (circles with error bars) by using Eq. (1231) . scatter around the experi- 
mental (thick dash-dotted lines) or theoretical (thick and thin lines) CE results. The results 
of calculations with the MCE-LNBCS (squares) and MCE-LNSCQRPA (triangles) by using 
the same definition ( 123]) and SS = 0.5 also describe well these values. The results for MCE 
entropies in Figs. [TJ and |3] show the importance of the effect beyond the quasiparticle mean 
field included in the self-consistent coupling QRPA vibrations. In fact, the MCE-LNSBCS 
results for the entropy clearly underestimate the experimental values. The discrepancy with 
the MCE-LNSCQRPA results increases with E* to reach about 20% at E* = 20 MeV. 



C. Level density 

The level densities obtained within the CE-LNSCQRPA using Eq. ([25]) and MCE- 
LNSCQRPA using Eq. (124"]) are plotted in Fig. H] as functions of excitation energy E* 



in comparison with the experimental data 
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2Q( Pobs{£) = Po x e w[Sobs(£)]- In. the latter 



Po is a normalization factor, which should be put equal to 1/SS according Eq. ( 1271) . How- 
ever, because of fluctuations in level spacings, which make the entropy sensitive to SS, the 
authors of Ref. [h], Q chose the values of po to obtain entropy S Q b s = at T = 0. In this 
way the value of p is set to 1.5 MeV" 1 for 94 - 98 Mo 2(J and 3 MeV" 1 for 162 Dy and 172 Yb 



19j. Figure H] shows that the level densities obtained within the MCE-LNSCQRPA offer the 
best fit to the experimental data for all nuclei under consideration. The results obtained 
within the CE-LNSCQRPA are closer to the experimental data for 94 ' 98 Mo at E* < 4 MeV, 
whereas at higher E* the MCE-LNSCQRPA offers a better performance. The S shape in 
the MCE-LNSCQRPA level density at low E* might have come from the fixed value of the 
energy interval SS, within which the levels are counted, according to the definition (T22]) . 
whereas the denominator in the definition of the CE level density [at the right-hand side of 
Eq. (|25p ] depends on E*. A larger value SS at E* < 4 MeV would eventually increase the 
MCE-LNSCQRPA level density, improving the agreement with the observed level density 
in this region, but there is no physical justification for doing so. The discrepancy between 
the CE-LNSCQRPA and experimental results seems to be larger and increases with E* for 
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FIG. 4. (Color online) Level densities as functions of E* obtained within the CE-LNSCQRPA (solid 
line) and MCE-LNSCQRPA (triangles) versus the experimental data (circles with error bars) for 
94 Mo (a), 98 Mo (b), 162 Dy (c), and 172 Yb (d). 

162 Dy and 172 Yb. This might be due to the absence of the contribution of higher multipolar- 
ities such as dipole, quadrupole etc., which are not included in the present study and may 
be important for rare-earth nuclei. On the other hand, the use of SCQRPA plus angular 
momentum [27(, discussed previously, may also improve the agreement. 



IV. CONCLUSIONS 



The present article applies the canonical and microcanonical ensembles of the LNBCS 



and LNSCQRPA approaches, derived in Ref. 22], to describe the thermodynamic properties 
as well as level densities of several nuclei, namely 94 ' 98 Mo, 162 Dy and 172 Yb. The results 
obtained show that the CE(MCE)-LNSCQRPA describe quite well the recent experimental 
level densities and the thermodynamic quantities extracted for these nuclei by the Oslo 
group 
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2l| . It confirms that the SN phase transition is smoothed out in nuclear systems 



due to the effects of quantal and thermal fluctuations leading to the nonvanishing pairing 
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gap at finite temperature even in heavy nuclei 3-8]. The discrepancy between the heat 
capacities obtained within the two different experimental works, which extrapolate the same 
experimental level density to different excitation energies, are also discussed. The heat 
capacities obtained within the CE-LNBCS(LNSCQRPA) for all nuclei show a pronounced 
oeak at T ~ T c , whereas the results extracted from the same experimental data by Refs. 



20| and 21] show different behaviors. The better agreement between the predictions of 



our approaches as well as those of the FTQTMC and the results of Ref. |2l| gives a strong 
indication to the fact that, to construct an adequate partition function for a good description 
of thermodynamic quantities, the measured level density should be extended up to very high 
excitation energy E* ~ 180 MeV or 200 MeV. The small differences between the CE(MCE)- 
LNBCS(LNSCQRPA) results and the experimental data might be due to the absence of the 
contribution of higher multipolarities such as dipole, quadrupole etc., which are not included 
in the present study. In order to tackle this issue, the LNSCQRPA plus angular momentum 
271 ] should be used and extended to included also the multipole residual interactions higher 
than the monopole pairing force. This task remains one of the subjects of our study in the 
future. 
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Appendix A: MCE results within the Richardson model 

The CE-LNBCS and CE-LNSCQRPA has been tested within the Richardson model in 
Ref. [22| and the results obtained are found in very good agreement with the exact solutions 
whenever the latter are available. In order to have more convincing evidences on the accuracy 
of present approaches, we show in Fig. [5] the MCE entropies and level densities obtained 
within the MCE-LNBCS and MCE-LNSCQRPA versus the exact ones for the Richardson 
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FIG. 5. (Color online) MCE entropies and level densities as functions of E* obtained within the 
MCE-LNBCS (squares), MCE-LNSCQRPA (triangles) versus the exact results for the Richardson 
model (circles) with N = Q = 14 and G = 1 MeV. Results obtained by using the energy bin 8£ = 
1 MeV are shown in (a) and (b), whereas those obtained by using 5£ = 5 MeV are shown in (c) 
and (d). Lines connecting the squares and triangles are drawn to guide the eye. 



model with N = £1 = 14 and G — 1 MeV. Two different values of energy interval 5£, namely 
S£ = 1 MeV (left panels) and SS = 5 MeV (right panels) are used in calculations. This figure 
shows that the MCE-LNSCQRPA always offers the best fit to the exact results, whereas 
the MCE-LNBCS underestimates the exact ones. The decreasing of the entropy as well as 
level density for the case with small value of 6£ = 1 MeV shown in Figs. |5] (a) and |5] (b) 
is due to the small configuration space with N = Q = 14 in the present case. This feature 
is ultimately related to the problem of using thermodynamics in very small system with 
discrete energy levels, where the temperature may decrease with increasing the excitation 



energy £* (See Fig. 2 of Ref. js]). This shortcoming can be effectively overcomed by using a 
larger 5£ = 5 MeV. As the result, the entropy and level density increase with increasing £* 
as shown in the right panels of Fig. [51 although there is no physical justification for using 
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such a large value of 58. 
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